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Abstract. These lecture notes are meant to provide a pedagogical introduction, and 
present the latest theoretical and experimental developments on the physics of vortices 
in type II superconductors. 



1 Why and what in vortices 

The discovery of high Tc superconductors has shattered the comforting sense 
of understanding that we had of the phase diagram and physical properties of 
type II superconductors, and in particular of the mixed (vortex) phase in such 
systems. Indeed, it was well known since Abrikosov [Q that above H c \ the mag- 
netic field penetrates under the form a vortex, made of a filament of radius £ 
(the coherence length) surrounded by supercurrent screening the external field 
running over a radius A (the penetration length). Because of the repulsion be- 
tween vortices due to supercurrents (see Fig. ^), the naive idea is that vortex 
will form a perfect triangular crystal (the Abrikosov lattice) . This has led to the 
phase diagram shown in Fig. [l], that has been the cornerstone of our understand- 
ing of all type II superconductors for more than three decades [^|,^| . However in 
high Tc, one could reach much higher temperatures, and it was soon apparent 
that some of the physics linked to the existence of the thermal fluctuations and 
disorder was overlooked. This led to a burst of investigations, both theoretical 
and experimental, to understand the physical properties of such vortex matter. 
Of course, high Tc were not the only field of investigations and low Tc supercon- 
ductors were reexamined as well, now that we knew what to look for in them. 

Indeed the vortex phase provides an excellent system for both the funda- 
mental researcher and one in search of physics with useful applications. From 
the fundamental point of view, vortex matter provides a unique system where 
one can study a crystal, in which one can vary the density (the lattice spacing) 
at the turn of a knob (simply by varying the magnetic field). In addition, be- 
cause this crystal is embedded in a "space" with a much finer lattice constant 
(the real atomic crystal), it can be submitted to various perturbations such as 
disorder, difficult to investigate in normal crystals. This provides thus a unique 
opportunity to study the combined effect of disorder and thermal fluctuations 
on a crystal. From a more practical point of view, if one has vortices in a su- 
perconductor and passes a current J in them, the current will act as a force on 
the magnetic tubes that are the vortices, and they will start to slide. The sliding 



2 T. Giamarchi and S. Bhattacharya 



V 



y 



«1 - 10000 A 

• 5 ~ioA 

•Hcl ~ 100 G 

* Hc2 - 300 T 



TTc2 



Hcl 



j i Abrikosov Vortex Crystal 




Normal 


p. > 




^J) (¥} O 




© -0 




Mcissncr 






— ► 



Tc T 



-8 Eh 



L- 



/Toe, 



Log(rl^) r<A 

e 



Fig. 1. The structure of a vortex, with a core size £ and supercurrent running over a 
radius A the penetration length. Typical values for high Tc are given. Due to super- 
currents vortex repel with a force F. This leads to the naive phase diagram where a 
crystal of vortex exists between H c \ and H C 2. 



in turn generates an electric field E, which means that the superconductor is 
not superconducting any more, due to the motion of vortices (see Fig. ^). The 




Fig. 2. A current J driven in a superconductor exert a force on the vortices, making 
them slide sideways with a velocity v. Since the vortices are flux tubes this generates 
an electric field E in the direction of the current and thus to a finite (and rather bad) 
resistance. 



resistance is rather poor and is simply p = p n {H j H C 2) (Bardeen-Stephen resis- 
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tance). In order to get a good superconductor, it is thus necessary to prevent the 
vortices from moving by pinning them. Hence the strong practical incentive to 
understand the properties (both static and dynamics) of vortices in the presence 
of disorder. 

A study of vortices prompts for several questions, that we will try to address 
in these notes: 

• What is the effect of disorder on the Abrikosov vortex crystal ? 

• How to describe the vortex phase ? Does one needs the full Ginzburg-Landau 
description or can one use a simplified description modelling vortices as 
elastic spaghettis ? 

• What is the phase diagram and the static physical properties of the vortices 
in presence of disorder and thermal fluctuations ? 

• What are the dynamical properties ? Is there a linear response and more 
generally what is the I — V characteristics ? 

• What is the nature of the vortex system when it is in motion ? 

• Are there links with other physical systems exhibiting the same competition 
between crystal order and disorder ? 

Of course, these questions have been examined intensively in the last 30 
years, and represent an impressive body of research work. So in these few pages 
we have to make a choice. Reviews already exist on vortices QJlJltO so we wiU 
try to present in these notes the basic ideas enabling the reader to understand 
the concept behind the variety of studies, and then bring the reader up to date 
with the recent theoretical and experimental developments that are left out of 
the previous reviews. The plan of the paper is as follows: in Section || we discuss 
an elastic description of the vortices, and the issue of lattice melting. In Section 
H we discuss the effects of disorder, introduce the basic lengthscales and physical 
concepts and discuss the limitations of the previously proposed solutions to 
tackle this problem. In Section ^ we discuss the recent theory of the Bragg glass 
and compare it with the host of experimental data on the statics of the vortex 
lattice. In Section |^ we discuss the more complicated issue of the dynamics. 
Finally conclusions, perspectives and contact with other physical systems can 
be found in Section |[ 

2 Elastic description of vortices 

A way to get a tractable description of the vortex lattice is to ignore the micro- 
scopic aspects of the superconducting state and the Ginzburg-Landau description 
of the vortex, and simply consider the vortices as an elastic object. The core is 
like a piece of string and the supercurrents provide the repulsive (elastic) forces. 
Of course such a description is a simplification and depending on the problem, 
it will be necessary to check that important physics has not been left out in the 
process. However, such a description has the advantage of being simple enough 
so that additional effects such as disorder can be included, and retains in fact 
most of the interesting physics for amacroscopic description of the vortex lattice 
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(phase diagram, imaging, transport). Another advantage is that this allows us 
to make contact with a large body of related problems as will be discussed in 
section ||. 

We thus describe the vortex system as objects having an equilibrium posi- 
tion (on a triangular lattice for the vortices, but this is of course general) 
and a displacement Ui compared to this equilibrium position. Ui is a vector with 
a certain number n of components. The vortices being lines n = 2, since dis- 
placements are on the plane perpendicular to the z axis. The elastic Hamiltonian 
is ^ 

H = ^^^ c ^{<l) u oL{q)up{-q) (1) 

a/3 q 

where a = x,y,z are the spatial coordinates. The c a p(q) are the elastic constants. 
The fact that they have a non trivial dependence on q comes from the long 
range nature of the forces between the vortices. The c can be computed from the 
microscopic forces between vortices. Standard elasticity corresponds to c(q) = 
cq 2 . Such a behavior will always be correct at large distance (small q) since the 
forces have a finite range A. In ([!]) various physical process have in principle to 
be distinguished and correspond to different elastic constants. This corresponds 
to bulk, shear and tilt deformations of the vortex lattice as shown on Fig. 
Although these different elastic constants can be widely different in magnitude 
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Fig. 3. Three different deformations of the vortex lattice (compression (bulk), shear 
and tilt) correspond to three different elastic constants, respectivelly called en, cg6 and 
C44. 



(for example bulk compression is usually much more expensive than shear), this 
is a simple practical complication that does not change the quadratic nature of 
the elastic Hamiltonian. Such a description is of course also valid for anisotropic 
superconductors (such as the layered High Tc ones), provided that the anisotropy 
is not too large. If the material is too layered then it is better to view the vortices 
as pancakes living in each plane and coupled by Josephson or electromagnetic 
coupling between the planes as shown in Fig. ^. For moderately anisotropic 
materials viewing the stack of pancakes as a vortex line is however enough. In 
this notes we will stick to this description. 

The melting of the vortex lattice can easily be extracted from the elastic 
description. Although a detailed theory of melting is still lacking, one can use a 
basic criterion, known as Lindemann criterion that states that the crystal melts 
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Fig. 4. In very anisotropic materials, it is best to see a vortex line as a stack of 
pancakes vortices coupled by electromagnetic or Josephson coupling. For most cases 
however a line description will be sufficient. 



when the thermally induced displacements a particle in the crystal becomes 
some sizeable fraction of the lattice spacing. On a more formal level the melting 
is given by 



(« > = 1 



2 „2 



(2) 



which defined the "effective" (thermal) size of the particle (also known as Linde- 
mann length) . The proportionality constant that reproduces correctly the melt- 
ing is empirically determined to be Cl ~ 0.1 — 0.2. A simple calculation based 
on the elastic the elastic description thus gives 
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(3) 



showing that the melting temperature goes down as the lattice spacing goes down 
(or the magnetic field up). Of course the full quantitative study for the vortex 
lattice should be done with the full fledged elastic Hamiltonian (including bulk, 
shear, tilt), but the main conclusions remain unchanged ||||. This leads to the 
first modification of the naive Abrikosov phase diagram taking into account the 
melting shown in Fig. ||. Close to H c \ the elastic constants drop down (since the 
vortices get separated by more than A the force between them is exponentially 
small) , and from (||) the crystal also melts, leading the reentrant behavior shown 
in Fig. | 

Early experimental studies p^ , |ll|,^2|Jl^| of the melting transition were based 
on transport measurements. The sudden onset of an ohmic resistance was argued 
to signify melting and its location in the (H,T) space is the locus of the melting 
phase boundary. Typical experimental results are shown in Fig. ^, where the 
onset is characterized by a pronounced knee in the resistance. With increasing 
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Tc T 

Fig. 5. Thermal fluctuations induce melting of the vortex lattice. On the left the full 
melting curve is shown. For high Tc, in practice one is often far from both H c \ and 
H C 2 leading to the apparent phase diagram shown on the right. 
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Fig. 6. T-dependence of ohmic resistance at various field values. The sudden onset 
marked by a pronounced knee marks the melting transition. The inset shows a plot 
of the temperature derivative of the resistance which, surprisingly mimics the specific 
heat jump across the transition Jl3j 
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external field, the onset moves to lower temperatures and eventually broadens 
considerably. This implies that at sufficiently large fields, the sharp first order 
transition crosses over to a more continuous second order like transition, an ef- 
fect expected to be the result of disorder (see later). At lower fields and higher 
temperatures, the melting transition is closely approximated by the disorder-free 
case, discussed above. Later experiments performed on thermodynamic quanti- 
ties such as the magnetization and specific heat confirmed the first order nature 
of the melting transition, at least for weak disorder. Fig. [?] shows an experi- 
mental phase diagram of BSCCO obtained by local magnetization using a novel 
hall-bar technique 0. Figure || shows typical measurements of the jump in the 
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Fig. 7. Melting phase diagram obtained from the jump in local induction for both 
isothermal and isofield data fh4| . Note the similarity of the phase boundary with the- 
oretical expectations in Fig. plabove. 




local induction across this transition, by either changing temperature at a fixed 
field or changing field at a fixed temperature. In both cases, a positive jump in 
B is observed in going from the solid to the liquid phase. Using the Clapeyron- 
equation : 

dH m = _AS_ 

dT AM [ ' 

where AS is the entropy change and AM is the change in magnetization (den- 
sity of vortices). The experimental phase diagram is consistent with theoretical 
expectations as shown in Fig. [5] above. The negative slope of the melting curve 
is consistent with an increase of density in the liquid phase, curiously akin to a 
"water-like" melting phenomenon. Very little experimental work is available on 
the low field reentrant branch of the melting phase boundary. At low fields the 
intervortex interaction is weak and the effects of disorder dominate. Reentrant 
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Fig. 8. Typical jumps in local induction for isothermal and isofield measurements from 
Hall bar method J14| . Both the solid and liquid phases have no measurable pinning in 
the bulk and the data are reversible. 



phenomena in peak effect (see Section 4.4) have been somewhat widely observed 
and is thought to be dominated by effects of disorder, rather than thermal fluc- 
tuations. 

Finally, besides the phase diagram, what are the physical quantities that one 
can in principle compute and that are directly connected to experiments ? The 
first important information is the relative displacements correlation function 



B(r) = <[«(r)-u(0)P) 



(5) 



where () is the average over thermal fluctuations and 7TT is the average over 
disorder (if need be) . (0) indicates how the displacements between two points in 
the system separated by a distance r grow (see Fig. ^J). Decoration experiments 




Fig. 9. The relative displacement correlation function B(r) measuring the displace- 
ments between two points separated by a distance r is directly measured in experiments 
such as decoration experiments. 



provide a direct measure of this correlation function as we will see. In a perfect 
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crystal B(r) = 0, whereas both thermal fluctuations and disorder will make the 
displacements grow. How B{r) grows tells us whether the system is well ordered 
or not. In a good crystal B(r) will saturate to a finite value whereas it will grow 
unboundedly if the perfect positional order of the crystal is destroyed. 

Another important quantity, directly measures in diffraction experiments is 
the structure factor 

S(q) = (p(q)p(-q)) (6) 

In a perfect crystal this consists of Bragg peaks at the vectors K of the reciprocal 
lattice of the crystal. If one considers one such peak its shape (as shown on 
Fig. |ni| ) is the Fourier transform of the positional correlation function 
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Fig. 10. The structure factor measured in diffraction experiments such as neutrons 
or X-rays. The shape of a peak is given by the Fourier transform of the positional 
correlation function C(r) (see text). 



C{r) = ( e iKu(r) e -iKu(^ ( 7 ) 

Thus in a perfect crystal C(r) = 1 and the Fourier transform is a 6(q) Bragg 
peak. If there are only thermal fluctuations C(r — > oo) — > Cste (in fact C(oo) = 
e -K i T /2y p our j er transform is still a S peak but with a reduced weight 
which is simply the Debye Waller factor. The faster C(r) decreases the more 
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disordered is the crystal. If C(r) decreases exponentially to zero with a char- 
acteristics lengthscale R a , the peak in the structure factor is some patatoidal 
(lorentzian like) shape with a width R^ 1 indicating that the perfect transla- 
tional order is lost. Although it is not always true (it is only true for gaussian 
fluctuations such as thermal fluctuations) a rule of thumb is 

C(r) ~ e -^M/ 2 (8) 

showing quite logically that the faster the relative displacements grow the more 
positional order (measured by the Bragg peaks) is destroyed in the system. 

3 Disorder, basic lengths and open questions 

The next task is to consider the effects of disorder. In real systems, disorder exists 
in all varieties in the underlying atomic crystal: vacancies, interstitials, lattice 
dislocations, grain boundaries, twin boundaries, second-phase precipitates, etc. 
In high quality single crystals it is possible to limit dominant disorder to point 
like impurities. Additionally, point like disorder can be, and sometimes is, in- 
tentionally added to systems in the form of electron irradiation in the case of 
the cuprates and/or substitutional atomic impurities of various kinds in low Tc 
systems. More artificial disorder can be introduced for example by heavy ions 
irradiation that produce columns of defects in the material. We will briefly men- 
tion the consequences of such artificial disorder in Section |^, but most of these 
notes will be devoted to the effects of point like impurities. 

Such impurities (as shown on Fig. Ill]) can be modeled by a random potential 




Fig. 11. Point like impurities (small circles) act as pinning centers for the vortices. 
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V(r) coupled directly to the density p(r) = ^ ■ 5(r — — m) of vortex lines 

H dis = / d d rV(r)p(r) (9) 



In principle one has "just" to add (||) to (|j) and solve. Unfortunately the coupling 
of the displacements to disorder is highly non linear (since it occurs inside a 
S function), and thus this is an horribly complicated problem. Physically this 
traduces the fact that there is a competition between the elastic forces that want 
the system perfectly ordered and the disorder that let the lines meander. This 
competition is bound to lead to complicated states where the system tries to 
compromise between these two opposite tendencies. 

In order to understand the basic physics of such problem a simple (but ground 
breaking !) scaling argument was put forward by Larkin fl5|| . To know whether 
the disorder is relevant and destroys the perfect crystalline order, let us assume 
that there exists a characteristic lengthscale R a for which the relative displace- 
ments are of the order of the lattice spacing u(R a )—u(0) ~ a. If the displacements 
vary of order a over the lengthscale R a the cost in elastic energy from (|]) is 

^-V (io) 

by simple scaling analysis. Thus in the absence of disorder minimizing the energy 
would lead to R a = oo and thus to a perfect crystal. In presence of the disorder 
the fact that displacements can adjust to take advantage of the pinning center 
on a volume of size R d allow to gain some energy. Since V is random the energy 
gained by adapting to the random potential is the square root of the potential 
over the volume R d , thus one gains an energy from (^|) 

- VR d J 2 Po (11) 

Thus minimizing ( |To| ) plus ( |TT[ ) shows that below four dimensions the disorder 
is always relevant and leads to a finite lengthscale 

/ cV N 1/(4 ~ d) 

a - 



Ra ~ a [v^ ) (12) 

at which the displacements are of order a. The conclusion is thus that even an 
arbitrarily weak disorder destroys the perfect positional order below four dimen- 
sions, and thus no disordered crystal can exist for d < 4. This is an astonishing 
result, which has been rediscovered in other context (for charge density waves 
R a is known as Fukuyama-Lee Jl6| length and for random field Ising model this 
is the Imry-Ma length fl7|| ). Of course it immediately prompt the question of 
what is the resulting phase of elastic system plus disorder ? 

Since solving the full problem is tough another important step was made by 
Larkin |l5|,^8). For small displacements he realized that (||) could be expanded 
in powers of u leading to the simpler disorder term 

fflarkin = / d d rf(r)u(r) (13) 
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where f(r) is some random force acting on the vortices. Because the coupling 
to disorder is now linear in the displacements the Larkin Hamiltonian is exactly 
solvable. Taking a local random force f(r)f(r') = AS(r—r') gives for the relative 
displacements correlation function and structure factor 

B(r) = SthcrmalM + 4^ (14) 

C(r) = e~ K ' B ^' 2 ~ e~ r " d (15) 

where -Bthormai are the displacements in the absence of disorder due to thermal 
fluctuations (which remain bounded in d > 2 and are thus negligeable at large 
distance compared to the disorder term) . Thus the solution of the Larkin model 
confirms the scaling analysis: (i) displacements do grow unboundedly (as a power 
law) and thus perfect positional of the crystal is lost; (ii) the lengthscale at which 
the displacements are of the order of the lattice is the similar to the one given by 
the scaling analysis. In addition the Larkin model tells us how fast the positional 
order is destroyed: the displacements grow as power law thus the positional is 
destroyed exponentially fast, leading to peaks in the structure factor of width 
i?" 1 . However all these conclusions should be taken with a grain of salt. Indeed 
the Larkin model is an expansion in powers of u, and thus cannot be valid at large 
distance (since the displacements grow unboundedly the expansion has to break 
down at some lengthscale). What is this characteristic lengthscale ? A naive 
expectation is that the Larkin model cease to be valid when the displacements 
are of order a i.e. at r = R a . In fact this is too naive as was noticed by Larkin 
and Ovchinikov. To understand why, in a transparent way, let us rewrite the 
density of vortices in a way more transparent than the original form 

p(x) = J2S(r-R° i -u i ) (16) 

i 

Taking the continuum limit for the displacements in ( |l6| ) should be done with 
care since one is interested in variations of the density at scales that can be 
smaller than the lattice spacing. A very useful way to rewrite the density is 
[[H|§ (see also @||): 

p(r) =p - Po V ■ u(r) + p Q £ e iK ^~ u ^ (17) 

K 

which is a decomposition of the density in Fourier harmonics determined by the 
periodicity of the underlying perfect crystal as shown on Fig. [l^. The sum over 
the reciprocal lattice vectors K obviously reproduces the <5 function peaks of the 
density (|l6|). If one considers particles with a given size (for vortices it is the 
core size £) then the maximum K vector in the sum should be 

K max ~ 2tt/( (18) 

in order to reproduce the broadening of order £ of the peaks in density. This 
immediately allows us to reproduce the Larkin model by expanding (0) using 
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Fig. 12. Various harmonics of the density. If one is only interesting in variations of 
the density at lengthscales large compared to the lattice spacing one has the standard 
"elastic" expression of the density in terms of the displacements. In the vortex system 
however on has to consider variations of density at lengthscales smaller than the lattice 
spacing and higher harmonics are needed (lEjJ2C|] . 



in powers of u. Clearly the expansion is valid as long as K max u <C 1 This 
will thus be valid up to a lengthscale R c such that u(R c ) is of the order of 
the size of the particles £. Note that this lengthscale is different (and quite 
generally smaller) than the lengthscale R a at which the displacements are of the 
order of the lattice spacing. The Larkin model cease to be valid way before the 
displacements become of the order of a and thus cannot be used to deduce the 
behavior of the positional order at large length scale. In addition it is easy to 
check that because the coupling to disorder is linear in the Larkin model, this 
model does not exhibit any pinning. Any addition to an external force leads to a 
sliding of the vortex lattice. It thus seems that this model is not containing the 
basic physics needed to describe the vortex lattice. In a masterstroke of physical 
intuition Larkin realized that the lengthscale at which this model breaks down is 
precisely the lengthscale at which pinning appears 111] . The lengthscale R c is thus 
the lengthscale above which various chunks of the vortex system are collectively 
pinned by the disorder. A simple scaling analysis on the energy gained when 
putting an external force 



allows to determine the critical force needed to unpin the lattice. Assuming that 
the critical force needed to unpin the lattice is when the energy gained by moving 
due to the external force is equal to the balance of elastic energy and disorder 




K 



(19) 




(20) 



|£ 2 i?J? 2 , one obtains 




(21) 
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This is the famous Larkin-Ovchinnikov relation which allows to relate a dynam- 
ical quantity (the critical current at T — needed to unpin the lattice) to purely 
static lengthscales, here the Larkin-Ovchinikov length at which the displace- 
ments are of the order of the size of the particle. Let us insist again that this 
lengthscale controling pinning is quite different from the one R a at which dis- 
placements are of the order of the lattice spacing at that controls the properties 
of the positional order. 

The lack of efficiency of the Larkin model to describe the behavior of displace- 
ments beyond R c still leaves us with the question of the nature of the positional 
order at large distances. However, extrapolating naively the Larkin model would 
give a power law growth of displacements. Such behavior is in agreement with 
exact solutions of interface problems in random environments (so called random 
manifold problems) and solutions in one spatial dimension. It was thus quite 
naturally assumed that an algebraic growth of displacements was the correct 
physical solution of the problem, and thus that the positional order would be 
destroyed exponentially beyond the length R a . This led to an image of the dis- 
ordered vortex lattice that consisted of a crystal "broken" into crystallites of 
size R a due to disorder. To reinforce this image (incorrect) "proofs" were given 
[ p3| to show that due to disorder dislocations would be generated at the length- 
scale R a (even at T — 0) further breaking the crystal apart and leaving no hope 
of keeping positional order beyond R a . A summary of this (incorrect) physical 
image is shown on Fig. |l3|. 




Fig. 13. The (incorrect) physical image that was the commonly accepted view of 
what a disordered elastic system would look like. The crystal would be broken into 
crystallites of size R a by the disorder. Dislocations would be generated by the disorder 
at the same lengthscale. 

If one believes that the positional order is lost and dislocations are sponta- 
neously generated one can wonder whether an elastic description of the vortex 
lattice is a good starting point. An intermediate attitude is to consider that such 
a description is useful at intermediate lengthscale and can be used to obtain the 
pinning properties (since they are controlled by lengthscales below R a ) or as a 
first step in absence of a better description [Q. A more radical view is to con- 
sider that since positional order would be lost at large length scale it is best to 
ignore it from the start and that an elastic description of the vortex lattice is a 
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bad starting point: it is much better to ignore positional order altogether and to 
focus on the phase of a vortex |^4|,^3| . The effect of disorder is thus introduced 
by a random gauge field destroying the phase coherence between the vortices. 
The system is then described by a random phase energy 

H = co <& ~ 4>j ~ Mj ) ( 22 ) 

ij 

With certain assumptions on the properties of the gauge field (essentially that 
A = oo) the idea is that the solid vortex phase will be transformed into a glassy 
Gauge glass (called the vortex glass), leading to the phase diagram shown in 
Fig. |l4j. The vortex glass phase has a continuous transition, with a divergent 




Fig. 14. In a vision where positional order is ignored, and vortices are described 
by they phases, the solid phase is transformed into a vortex glass phase, having a 
continuous transition (with scaling) towards the liquid phase [ p4p3| . 



lengthscale, towards the liquid. It thus exhibit scaling at the transition. It was 
also suggested that inside the glass phase there should be no linear response to 
an applied current [H . We will come back to this point in Section [| 

Although this description of the vortex lattice/ vortex glass phase was very 
successful in the beginning, it started to run into serious problems as both the 
experiments and the theory were refining. Among the experimental problems 
one would notice (the corresponding data will be presented in the next sections) 

• The transition between the solid (vortex glass ?) and the liquid was shown 
to be discontinuous by various measurements. Specific heat measurements 
have now proved that this transition is first order. 

• Decoration experiments were seeing very large regions free of dislocations 
and showing a very good degree of positional order. This did not seem to fit 
well with the idea that disorder would strongly affect the positional order. 

• Neutron scattering was exhibiting quite good Bragg peaks, again showing 
stronger positional order than naively anticipated. 

• The phase diagram seemed more complicated than the one shown in Fig. |l4]. 

On the theoretical front two main points were raised: (i) the gauge glass 
model was shown to have no glass transition in d = 3 for realistic values of the 
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parameters (such as a finite A) p5[ . (ii) The important suggestion was made, us- 
ing scaling arguments [ p2[]2f| and then firmly established in detailed calculations 
[ p7p9| that provided dislocations were ignored displacements in vortex lattices 
were growing much more slowly than a power law (logarithmically) . 

These experimental facts and theoretical points suggested that the effect of 
disorder on the vortex lattice could be less destructive than naively anticipated. 
They thus strongly prompted for an understanding of the physical properties, 
such as the positional order, stemming from the elastic description. They also 
made it mandatory to resolve the issue of the asserted (2(||23| ever presence 
of disorder induced dislocations, which would invalidate the elastic results and 
always destroy the positional order above R a . 



4 Statics : Experimental facts and Bragg glass theory 

To get a quantitative theory of the disordered system, and go beyond the simple 
scaling analysis, it is necessary to solve the full ([!]) plus (||) . Fortunately the the- 
oretical "technology" had developed tools allowing to obtain a rather complete 
solution of this problem [[l9|]2C|| . We describe the solution here and examine the 
consequences for experimental systems. 



4.1 Bragg glass 

The problem one needs to solve is (using the decomposition of density (|l 



R =\\ d d r(Vu) 2 + Po J2 f d d re lK{r - u{r)) V{r) (23) 
J K J 

Although we have written here the simplified form of the elastic hamiltonian the 
full one has to be considered but this does not change the method. One then 
gets rid of the disorder using replicas. After averaging over disorder the problem 
to solve becomes 

H = \ E / A (V« a ) 2 -DplY, E / ddr «*(tf(u a (r) - u b (r))) (24) 



where V(r)V(r') = DS(r — r'). One has thus traded a disordered problem for a 
problem of n interacting fields. The limit n — ► has to be taken at the end for 
the two problems to be identical. So far the mapping is exact, but (53) is still 
too complicated to be solved exactly. Two methods are available to tackle it: 
(i) a variational method; (ii) a renormalization group method around the upper 
critical dimension d = 4 (a 4 — e expansion). The renormalization method is 
relatively involved and we refer the reader to the various reviews and to ]20| , p8[ 
for more details and discussions. The variational method is simple in principle 
[ p9| , and has the advantage to give the essential physics. One looks for the best 
quadratic Hamiltonian 



E / d d qG-J{q)u a {q)u b {-q) (25) 
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that approximate (p4|). Hq leads to a variational free energy 

F var = F + (H - H ) Ho (26) 

that has to be minimized with respect to the variational parameters. The un- 
known Green's function G a b(q) are thus determined by 

r)F 

var =0 (27) 



dGab{q) 



This is nothing but the well known self consistent harmonic approximation. 
The technical complication here consists in taking the limit n — > [[30). The 
best variational parameters are the ones that break the replica symmetry, in a 
similar way than what happens in spin glasses. This is very comforting since we 
expect on physical grounds that the competition between the elasticity and the 
disorder causes a strong competition where the system has to find its ground 
state. It is thus quite natural that such a competition leads to glassy properties. 
This is what the solution of the problem confirms. A similar effect appears in the 
renormalization solution where a non-analyticity appears, signaling again glassy 
properties. The two methods thus agree quite well (not only qualitatively but 
also quantitatively). 

Let us now describe the full solution given in |l^,^| . One finds for the relative 
displacements correlation function the one shown on Fig. Hq. Three regime can 
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Fig. 15. The relative displacement correlation function B(r) as a function of distance. 
The lengthscales R c and R a define three regimes. Below R c one recovers the Larkin 
behavior. Between R c and R a the displacements still grow algebraically albeit with a 
modified exponent. Above R a the growth becomes logarithmic. 



be distinguished, separated by the characteristic lengthscales R c and R a . Below 
R c , one is in the Larkin regime, and the displacements grow as r' 4- ^/ 2 . Then 
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between R c and R a (i.e. when the displacements are between £ and a), each line 
wanders around its equilibrium position and the problem is very much like the 
one of a single line in a disordered environment, i.e. a random manifold problem. 
The growth of displacements is still algebraic, albeit with a different exponent 
v. Above R a however the displacements grow much more slowly and B(r) = 
-i4(jlog(r). The physics is easy to understand: because of the periodic nature of 
the system, each line can take care of the disorder around its equilibrium position. 
There is thus no interest for one line to make displacements much larger than 
the lattice spacing to pass through a particularly favorable region of disorder 
(this would be the case for a single line). For many lines it is better to let the 
line closeby take care of this region. This ensures that the total energy of the 
whole system is the lowest. As a results the displacements do not need to grow 
much above a. The variational method and RG techniques allow to compute 
the prefactor Ad and to address the issue of the positional order. The positional 
correlation function is simply given by fl9| , pof 

C7(r)~(iy (28) 

where r\ is an exponent independent of the strength of disorder or temperature 
(for example ?7 V ar = 1 in rf = 3). The physics described above, and obtained 
from the variational ansatz is totally generic. The alternate RG approach in 
4 — e does indeed recover identical physics, as was shown first in the Larkin 



and asymptotic regimes 19 20fl and more recently for the full B(r) [|28|. The 
variational approach even gives quite accurate values of the exponents themselves 
as can be checked by comparing with the RG. The quite striking consequence 
is that far from having the positional order destroyed by the disorder in an 
exponential fashion, a quasi-long range order (algebraic decay of positional order) 
exists in the system. Algebraically divergent Bragg peaks still exist as shown on 
Fig. [l6[ It is to be noted that this phase, which is thus practically as ordered 
as a perfect solid, is a glass when one looks at its dynamical properties. This 
is seen in the analysis by the presence of replica symmetry breaking in the 
variational approach or the existence of non-analyticities in the renormalization 
solution. From a physical point of view this means that the system has many 
mctastable states separated from its ground state by divergent barriers. As a 
result it exhibits pinning and non-linear dynamics (creep motion), as we will 
discuss in more details in Section || 

What about the argument that dislocations should always be generated at 
lengthscale R a ? In fact it was shown J20| that this argument which forgets 
the fact that the coupling of the displacements to the disorder is non linear 
is simply incorrect and that in fact in d > 2 dislocations are not generated by 
disorder provided the disorder is moderate (i.e. R a large enough compared to a). 
This means that the results given above corresponds to the true thermodynamic 
ground state of the system. Thus there exists a thermodynamically stable glassy 
phase, the Bragg glass, with quasi long range order (algebraically divergent Bragg 
peaks), perfect topological order (absence of defects such as dislocations). This 
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Fig. 16. The disordered system still has quasi-long range positional order and alge- 
braically divergent Bragg peaks. It is thus nearly as ordered as a perfect solid. Such 
a phase, which is a glass when one looks at its dynamics properties, and in addition 
possesses perfect topological order (no defects such as dislocations) is the Bragg glass 
phase po| . 



is a surprising result since one naively associate a glass with a very scrambled 
system. The Bragg glass shows that this is not the case and that one has to 
distinguish the positional properties from the energy landscape (or the dynamical 
ones). 

The existence of this Bragg glass phase has of course many consequences for 
the vortex systems, consequences that we now investigate and compare with the 
available experimental data on vortices. 



4.2 Positional order: Decorations and Neutrons 

The first consequence is the existence of the perfect topological order and the 
algebraic Bragg peaks. 

On the theoretical side, after the initial proposal [E0| the existence of al- 
gegraic bragg peaks and the absence of dislocations have been confirmed by 



further analytical results |3l| , |32| , |33| and numerical simulations []34 35 1. On the 
experimental side, direct structural information has been obtained from both 
real space studies such as magnetic decoration, scanning tunneling microscopy, 
Lorentz microscopy and reciprocal space studies by neutron diffraction, sum- 
marized in Fig. 0. The upper two panels show a decoration micrograph of a 
fairly ordered lattice in NbSe2 with a few defects |3(| and an STM micrograph 
| |37| of the same material with no defects. The Lorentz micrograph in a thin 
film of Nb |j8| produces a nearly amorphous vortex assembly while the neutron 
diffraction data on a single crystal sample of Nb |39| show Bragg peaks up to 
third order reflections, suggesting a very high degree of order. These results show 
that defect-free phases with Bragg reflections are experimentally observed, while 
highly defective or even amorphous or liquid like phases also exist. The task is 
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to find which parts of the phase space are occupied by each and what controls 
the phase transformations among them. 

The fourier transforms of the real space data [ p^|36|]4l| ] show very large re- 
gions free of dislocations and yield Bragg peaks routinely, suggesting a much 
stronger solid like order than would be expected naively from a vortex glass 
model. The situation with neutron diffraction is similar as shown on Figure [171 




Fig. 17. Vortex phase structure using various methods. The upper two panels show a 
decoration data on the left]36| and an STM data on the right in NbSe2j3^|. The lower 
two panels show a Lorentz holography micrograph of a thin film Nb|$8| and a neutron 
diffraction picture of a single crystal Nb[B9| 



Recently neutron data have provided a direct evidence of the presence of the 
Bragg glass phase. Indeed the power-law Bragg peaks shown in Fig. |l^ gives 
when convolved with the experimental resolution the result shown in Fig. fig . 
The width at half width of the observed peak is constant and determined by the 
experimental resolution, and the height is fixed by the position correlation length 
R a - Thus if disorder (or magnetic field - see in the next section) is increased 
the Bragg glass predicts that observed neutron peaks should collapse without 
broadening. This behavior, has been quantitatively tested on the compound 
(K,Ba)BiC>3 (see Fig. |l^) which has the advantage of being totally isotropic 



Vortex phases 21 



and thus avoid all complications associated with anisotropy such as possible 2D- 
3D crossovers. Peaks are seen to collapse without any broadening thus providing 
a direct evidence of the Bragg glass phase and its algebraic positional order. 



1.2 




-l/R„ -1/4 l/£ i/R a Q' =K <*° 



Fig. 18. (left) Bragg glass predictions for the angular dependence of the neutron 
diffracted intensity when a finite experimental resolution is taken into account. The 
arrows indicate the values of S(q) for q = and 1/R a respectively. If the disorder 
increases (Ra decreases), the height of the peak decreases as R2 but the peak does 
not broaden since the half width at half maximum is always given by the experimental 
resolution l/£ and not 1/R a - The height of the peaks gives a direct measure of the 
characteristic length R a . (right) Angular dependence of neutron intensity diffracted by 
a (K,Ba)Bi03 crystal p2| at the indicated applied fields. Those data show that the 
diffracted intensity (rocking curves) collapse without any broadening above 0.7T. 



4.3 Unified phase diagram 

Another striking consequences of the existence of the Bragg glass, is that it im- 
poses a generic phase diagram for all type II superconductors [ p0[ . Indeed the 
Bragg glass has no dislocations, thus if either thermal fluctuations or disorder 
are increased the Bragg glass should "melt" to a phase that contains defects . If 
thermal fluctuations increase this is the standard melting and leads to the liquid 
phase. Because the Bragg glass is nearly as ordered as a perfect solid one can 
expect it to melt though a first order phase transition. More importantly this 
"melting" of the Bragg glass can also occur because the disorder is increased 
in the system. For vortices increasing the field has a similar effect. Indeed the 
effective disorder in @ is Vpo, thus increasing the average density makes the 
disorder term stronger compared to the elastic term ([j]). Indeed for moderate 
fields the change in elastic constants due to the field is quite small. Thus in- 
creasing the field is like increasing disorder. One should thus have a "melting" 
transition of the Bragg glass (induced by the disorder) as a function of the field. 
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Close to H c i because the change of elastic constants is then drastic, one expects 
a similar transition. This leads to the universal phase diagram shown in Fig. [19. 



H,A 




Fig. 19. Universal phase diagram for type II superconductors. The Bragg glass "melts" 
due to thermal fluctuations of disorder induced fluctuations. This lead to transitions 
as a function of magnetic field. The left diagram is far from H c \ (more adapted to 
High Tc) whereas the insert shown the full diagram. The melting towards the liquid is 
expected to be first order. 



We first focus on the experimental results of the loss of Bragg glass order 
in the cuprate systems. A remarkable experimental determination of the phase 
diagram of BSCCO is shown in Fig. ^0|. Three phases are clearly identified in it : 
a quasi-ordered Bragg glass phase which "melts" by thermal fluctuations into a 
liquid and also "amorphizes" or "melts" by quenched disorder into a disordered 
solid. First we focus on structural evidence of the two phase transitions from 
neutron diffraction [H. In the Bragg glass phase one clearly observes the reso- 
lution limited Bragg peaks. Upon increasing temperature across the thermally 
driven transition the Bragg peaks lose intensity and become unobservable in the 
putative liquid phase. A similar loss of Bragg peak intensity is observed when 
the magnetic field is increased at fixed T across the Bragg glass to the disordered 
or vortex glass phase. Fig. |l| and Fig. ^2] show the results of neutron diffraction 
of BSCCO as the two phase boundaries are crossed, one from the ordered Bragg 
glass to the liquid and the other from Bragg glass to the disordered solid phase 
[ p4| . In both cases the Bragg reflections lose intensity and disappear at the phase 
boundary, in a way similar to the one discussed in Fig. 18. Due to limited neutron 
intensity, detailed studies of the disordered phase has not been performed in the 
cuprates. But recent neutron diffraction studies of low Tc system Nb f4j|,|39| with 
a short penetration depth have directly shown a transformation of bragg peaks 
to a ring of scattering, implying liquid like (amorphous) order in the disordered 
phase. The experimental phase behavior is thus entirely compatible with theoret- 
ical expectations in Fig. The comparison with the theoretical phase diagram 
identifies the quasi-ordered solid phase with the Bragg glass phase. The phase 
with dislocations is expected to correspond to the disordered solid phase. For 
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Fig. 20. Experimental phase diagram for the cuprate system BSCCO fljLEIi ]. At the 
thermally induced melting transition one sees thermodynamic anomalies as in Fig. |^. 
The disorder-driven transition is accompanied by the second peak or peak effect, shown 
later in Fig. 

BSCCO the position of the field melting line has been computed by Lindcmann 
argument or similar cage arguments p6p7| , p8| , p9{ and the value of the "melting" 
field is in good agreement with the observed experimental value. The distinction 
between the disordered solid phase and the liquid phase remains an experimen- 
tally open question for different systems with different types of disorder. For the 
very anisotropic BSCCO system there is also the question of the existence of 
additional phases. Structural results of the same kind are not available for the 
other common cuprate system, namely YBCO. However, thermodynamic data 
on the magnetization jump and entropy jump were measured. A composite of 
the data is shown in Fig. ^3] (see also ^o|) which demonstrates close agreement 
within the Clapeyron equation confirming the first order nature of the thermally 
driven melting transition. 

4.4 Second peak and peak effect 

Now we focus on the disorder driven phase transition shown above. The mag- 
netization hysteresis loop across this transition manifests a typical second peak 
effect, often called the fishtail effect due to its shape, as shown in Fig. |i| for 
BSCCO [ (43| . The sharp jump is magnetization at B sp marks a sudden increase 
in the irreversibility, i.e., a jump in critical current J c . In a naive view of the 




Fig. 22. Loss of Bragg intensity upon increasing the magnetic field at a fixed temper- 
ature 
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Fig. 23. Experimental data pl|j52] ] showing thermodynamic measurements from mag 
netization and calorimetry of the first order transition in YBCO 



Larkin scenario, this marks a sudden decrease in the correlation volume, i.e., a 
sharp loss of order, consistent with the neutron diffraction data shown above. 
Similar results were obtained for YBCO also, yielding a qualitatively similar 
phase diagram p3|,|54 55 . 



Peak effects are ubiquitous in low Tc systems and have been known for nearly 
four decades and provided the primary motivation for the Larkin- Ovchinnikov 
scenario of collective pinning. In these systems the peak effect |5(| usually occurs 
very close to the normal phase boundary, unlike in the cuprate systems where 
the fishtail anomaly occurs very far from it. Recent resurgence of activity on 
the peak effect phenomenon in low Tc systems also provide a phase diagram of 
ordered and disordered vortex phases not dissimilar to the cuprates. Due to the 
smallness of the thermal fluctuation effects, however, the peak effect transition 
often occurs in close proximity to the melting transition, or even coincident 
with it. Separating the two effects remains a matter of considerable controversy 
in the low Tc systems. In the popular low Tc system NbSe2 the peak effect 
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Fig. 24. Typical second magnetization peak data from the magnetic hysteresis loop 
in BSCCOp3|, the sudden enhancement of the diamagnetic signal at B sp marks the 
enhancement of critical current across the Bragg glass to vortex glass phase transition. 



phenomenon has been studied extensively in recent years [ p7[ . Fig. ^5] shows a 
typical set of data for the resistive detection of the peak effect where a rapid 
drop in the resistance at the peak effect boundary marks a sudden increase in 
the critical current analogous to the magnetization jump shown above. Direct 
structural evidence also clearly shows the amorphization of the Bragg glass phase 
with six fold symmetric Bragg spots to a ring of scattering at the peak in the 
elemental system Nb In NbSe2 the same peak effect is accompanied by a 
sharp change in the line shape as seen in the asymmetry parameter pd^ . These 
results are shown in Fig. |2(]. 

Of special importance is the clear experimental observation of a reentrant 
phase behavior for the NbSe2 system |^|5l). From the Meissner phase an in- 
crease in field shows two anomalies, first from a disordered phase into an ordered 
phase and then a reentry into a disordered phase just below the upper critical 
field. Especially significant is the pronounced shift in the order-disorder phase 
boundary with varying quenched disorder (addition of magnetic dopants). Fig- 
ure |27] shows the progressive shrinkage of the Bragg glass phase from both high 
field side as well as from the low field side as disorder is increased from sample A 
through sample C. These results are in excellent agreement with the theoretical 
discussions above. Direct structural evidence of amorphization in this system 
was obtained through muon spin relaxation experiments |59 that are entirely 
analogous to the results [[32| in the cuprate systems. 

Several questions remain open for the peak effect from the phase behavior 
point of view. In addition to the second magnetization peaks, a peak effect 
is often observed in YBCO very close to or even coincident with the melting 
transition [p3f . This suggests that the disordered solid phase may protrude as 
a sliver all around the Bragg glass phase f64[ or that there are two types of 
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Fig. 25. Typical observation of peak effect in transport ]5q |. The upper two panels 
show the transition in isothermal and isofield measurements. The lower left panel shows 
the independence of the locus on driving current suggesting a thermodynamic origin of 
the anomaly. The right panel shows the locus of the peak region bounded by the onset 
at H p i and the peak at H p . The close proximity to the upper critical field is typical of 
low T c materials.. 



peak effects, one associated with disorder induced melting and the other with 
the thermally induced melting transition. In what follows, we indeed show that 
the melting of the Bragg glass provides a natural explanation for the peak effect. 

How the melting of the Bragg glass signals itself ? To understand it let us look 
at the V — I characteristics, shown on Fig. |28|. The Bragg glass is collectively 
pinned so it has a small critical current J\ but very high barriers leading to 
practically no motion (hence no V) in the pinned phase (below Ji). On the 
other hand in the high field phase (with dislocations) or the liquid it is more 
easy to pin small parts leading to higher critical currents, but the pinning is not 
collective hence a much more linear response below J2 . The V — I characteristics 
thus cross at the melting of the Bragg glass This crossing leads for an 
apparent increase of the critical current close to the melting and thus to a peak 
effect in the transport measurements or to a second peak in the magnetization 
measurements . 



Recent simulations and reexamination of older experimental data |65 58 



are 



in excellent agreement with this I ~V crossing scenario at the peak effect near 
melting or the second magnetization peak. 
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Fig. 26. Experimental data on NhQ on the left and NbSe2 (59) on the right. The 
upper panel shows the peak effect from ac susceptibity for both. The lower left panel 
shows the structure from neutron diffraction in Nb showing the Bragg spots disappear 
into a ring of scattering. Muon spin relaxation data in NbSe2 on the right shows a 
sudden drop in the skewness parameter at the peak effect that marks the first order 
structural disordering of the lattice. 



5 Dynamics of vortices 

The competition between disorder and elasticity manifests also in the dynamics 
of such systems, and if any in a more dramatic manner. When looking at the 
dynamics, many questions arise. Some of them can be easily asked (but not 
easily answered) when looking at the V — I characteristics shown in Fig. ^9|. The 
simplest question is prompted by the T = behavior. Since the system is pinned 
the velocity is zero below a certain critical force F c . For F > F c the system moves. 
What is F c ? We saw that the scaling theory of Larkin and Ovchinikov relates it 
directly to the static characteristic length R c . Can one extract this critical force 
directly from the solution of the equation of motion of the vortex lines 

duAt) 5H 

This equation is written for overdamped dynamics, but can include inertia as 
well, rj is the friction coefficient taking into account the dissipation in the cores, 
F the externally applied force, and C a thermal noise. Can one obtain from this 
equation the velocity above F c ? The v — F curve at T = is reminiscent of the 
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Fig. 27. Evolution of the phase diagram [ pl| with increased disorder from sample A 
to C. As disorder increases, the collectively pinned lattice, i.e., the Bragg glass phase, 
shrinks from both above and below. 
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Fig. 28. The V — I characteristics for the Bragg glass (full line) and the high field 
phase or the liquid (dashed line) p?| ]. 
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Fig. 29. The V — I characteristics for a vortex system. This is in fact the velocity of 
the system v, in response to an external force F. 



one of an order parameter in a second order phase transition. Here the system 
is out of equilibrium so no direct analogy is possible but this suggests that one 
could expect v ~ (F — F C )P with a dynamical critical exponent (3. We will not 
investigate these questions here, because of lack of space and refer the reader 
to the above mentioned reviews and to |6(J for an up to date discussion of this 
issue and additional references. 

The second question comes from the T ^ curve. Well below threshold 
F<F C the system is still expected to move through thermal activation. What 
is the nature of this motion and what is the velocity ? If the system has a glassy 
nature one expects it to manifest strongly in this regime since thermal activation 
will have to overcome barriers between various states. We address this question 



in Section 5.1 
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Finally, there are many questions beyond the simple knowledge of the average 
velocity. One of the most interesting is the nature of the moving phase. If one is 
in the moving frame where the system looks motionless, how much this moving 
system resembles or not the static one. This concerns both the positional order 
properties and the fluctuations in velocity such as the ones measured in noise 



experiments. This is discussed in Section 5.2 



5.1 Creep 

Let us first examine the response of the vortex system to a very small external 
force. For usual systems we expect the response to be linear (leading to a finite 
resistivity for the system) . Indeed earlier theories of such a motion found a linear 
response. The idea is to consider that a blob of pinned material has to move in 
an energy landscape with barriers A as shown in Fig. BG. The external force 



A 




x 



Fig. 30. In the Thermally Assisted Flux Flow (TAFF) (67) a region of pinned material 
is considered as a particle moving in an energy landscape with barriers. This leads to 
an exponentionally small but linear response. 



F tilts the energy landscape making forward motion possible. The barriers are 
overcome by thermal activation (hence the name: Thermally Assisted Flux Flow) 
with an Arrhenius law. If the minima are separated by a distance a the velocity 
is 

v oc e ~^ A ~ Fa ^ - e -^ A+Fa ^ ~ e~ f3A F (30) 

The response is thus linear, but exponentially small. One thus recovers that 
pinning drastically improves the transport qualities of superconductors. For old 
superconductors (3 was small enough so that this formula was not seriously chal- 
lenged. However with high T c one could reach values of f} where it was clear that 
the TAFF formula was grossly overestimating the motion of vortex lines. 

The reasons is easy to understand if one remembers that the static system is 
in a vitreous state. In such states a "typical" barrier A does not exist, but barriers 
are expected to diverge as one gets closer to the ground state of the system. The 
TAFF formula is thus valid in system where the glassy aspect is killed. This is 
the case in the liquid where various parts of the system are pinned individually. 
When the glassy nature of the system persists up to arbitrarily large lengthscales 
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the theory should be accommodated to take into account the divergent barriers 
[ p4[|23[ |. This can be done quantitatively within the framework of the elastic 
description. In fact such a theory was developed before for interfaces |^,^| and 
then adapted for periodic systems such as the vortex lattice J7(],^2||. The basic 
idea is beautifully simple. It rests on two quite strong but reasonable assumptions 
: (i) the motion is so slow that one can consider at each stage the lattice as 
motionless and use the static description; (ii) the scaling for barriers which is 
quite difficult to determine is the same than the scaling of the minimum of energy 
(metastable states) that can be extracted again from the static calculation. If 
the displacements scale as u ~ L v then the energy of the metastable states (see 
(0)) scale as 

E ~ L d - 2+2v (31) 

on the other hand the energy gained from the external force over a motion on a 
distance u is 

E F = J d d xFu(x) ~ FL d+l/ (32) 

Thus in order to make the motion to the next metastable state one needs to 
move a piece of the pinned system of size 



(33) 



The size of the minimal nucleus able to move thus grows as the force decrease. 
Since the barriers to overcome grow with the size of the object the minimum 
barrier to overcome [assuming that the scaling of the barriers is also given by 
©) 

d-2+2v 

A(F)~(±) 2 ~" (34) 



leading to a velocity 

d-2+2v 

t)oce~ /3 (*) 2 ^ (35) 

This is a remarkable equation. It relates a dynamical property to static expo- 
nents, and shows clearly the glassy nature of the system. The corresponding 
motion has been called creep since it is a sub-linear response. Of course the 
derivation given here is phcnomenological, but it was recently possible to di- 
rectly derive the creep formula from the equation of motion of the system j7l|j6G|| . 
This proved the two underlying assumptions behind the creep formula and in 
particular that the scaling of the barriers and metastable states is similar. More 
importantly this derivation shows that although the formula for the velocity 
given by the phenomenological derivation is correct, the actual motion is more 
complicated than the naive phenomenological picture would suggest. Indeed the 
phenomenological image is that a nucleus of size L m i n moves through thermal 
activation over a length £^ lin , and then the process starts again in another part 
of the system. In the full solution, one finds that the motion of this small nucleus 
triggers an avalanche in the system over a much larger lengthscale [[36l . Checking 
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for this two scales process in simulations or actual experiments if of course a very 
challenging problem. 

The creep formula is quite general and will hold for interfaces as well as peri- 
odic systems. For periodic systems, dislocations might kill the collective behavior 
by providing an upper cutoff to the size of the system that behaves collectively 
(as if the system was torn into pieces). In the Bragg glass the situation is clear. 
Since there are no dislocations the creep behavior persists to arbitrarily large 
lengthscales. Since v — in the Bragg glass the creep exponent in ([35]) is fi = 0.5 
[22|2(|. What becomes of the creep when dislocations are present is still an open 
[ 72 and very challenging question. What is sure is that one can expect a weak- 
ening of the growth of the barriers or even their saturation, when going from the 
Bragg glass phase to the "melted" phase jlTj . This is at the root of the crossing 
of the I — V shown in Fig. Experimental verification of the creep effects 
postulated above have proved difficult due to the functional form of ( |35| ) where 
the power law appears in the exponentiated factor and requires data spanning 
many decades in the drive to determine the exponents with adequate precision. 
Transport experiments ]73| as well as magnetic relaxation experiments |74| have 
reported creep exponents compatible with the Bragg glass prediction, as well as 
weakening of the barrier growth when going to the disordered phase. But this is 
clearly a very challenging and difficult issue that would need more investigations. 



5.2 Dynamical phase diagram 

Let us now turn to the problem of the nature of the moving phase. 

This problem was directly prompted by experimental observations. Indeed 
early measurements of the dynamics showed dramatic effects near the peak ef- 
fect (see below), led to the construction of an experimental dynamical "phase 
diagram" |37| for the moving phases shown in Fig. |3l]. In addition dramatic 
evidence of the evolution of vortex correlations with driving force was obtained 
many years ago by neutron diffraction studies |75| . The Bragg peak in the pinned 
phase broadened significantly at the onset of motion showing a loss of order (or 
appearance of defects) and a subsequent healing of the Bragg peak at large 
drives showing a reentry into an ordered moving phase. It was thus necessary to 
understand the nature of the "phases" once the lattice was set into motion. 

One regime in which one could think to tackle this problem is when the lat- 
tice is moving at large velocities. Indeed in that case it is possible to make a large 
velocity expansion. Such expansion was performed with success to compute the 
corrections to the velocity due to pinning ]7(|(77]|, and get an estimate of the 
critical current. An important step was to use the large v method to compute 
the displacements [(78| . It was found that due to the motion the system averages 
over the disorder. As a result the system do not feel the disorder any more above 
a certain lengthscale and recrystalizes. The memory of the disorder would simply 
be kept in a shift of the effective temperature seen by this perfect crystal. This 
picture was consistent with what was shown to be true for interfaces, even close 
to threshold |79| (as shown on Fig. |3^) and thus provided a nice explanation for 
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Fig. 31. A typical dynamic phase diagram in NbSe2. Below the transition region a 
direct transition into a moving ordered phase is seen while in the peak regime a pro- 
nounced intermediate plastic flow regime is seen. At large drives above F cr , the moving 
ordered phase is established. In the equilibrium pinned liquid phase, the crossover cur- 
rent is immeasurably large. 

the recrystallization observed at sufficiently large velocities. Perturbation ap- 
proach along those lines has been extended in ^0|. However the peculiarities of 




Fig. 32. When an interface is in motion, it averages over the disorder. As a result the 
interface does not feel the disorder above a certain lengthscale R v and becomes flat 
again. 

the periodic structure manifests themselves again, and they does not follow this 
simple scenario, the way that the interfaces would. The crucial ingredient present 
in periodic system is the existence of a periodicity transverse to the direction of 
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motion. Because of this, the transverse components of the displacements still feel 
a disorder that is non averaged by the motion. As a result the large-w expansion 
always breaks down, even at large velocities and cannot be used to determine the 
nature of the moving phases. To describe such moving phase the most important 
components are the components transverse to the direction of motion (this is 
schematized in Fig. p3|). The motion of these components can be described by 
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Fig. 33. When in motion a periodic structure averages over the Fourier components of 
disorder along the direction of motion whereas the component perpendicular to motion 
remains unaffected. As a result the components of the displacements transverse to the 
direction of motion are the essential ingredients to describe a moving system. 



a quite generic equation of motion, as explained in [pl|j82| . The transverse com- 
ponents still experience a static disorder, and as a result the system in motion 
remains a glass (moving glass). 

In the moving glass the motion of the particles occurs through elastic chan- 
nels as shown in Fig. [54[ The channels are the best compromise between the 
elastic forces and the static disorder still experienced by the moving system. 
Like lines submitted to a static disorder the channels themselves are rough and 
can meander arbitrarily far from a straight line (displacements grow unbound- 
edly) . However although the channels themselves are rough, the particles of the 
system are bound to follow these channels and thus follow exactly the same 
trajectory when in motion. Needless to say the moving system (moving glass) 
is thus very different than a simple solid with a modified temperature where 
the particles would just follow straight line trajectories (with a finite thermal 
broadening) . 

When does this picture breaks down ? Clearly this should be the result 
of defects appearing in the structure. For example close to depinning, it was 
shown experimentally |57j (see Fig. |3l]) that some of the regions of the system 
can remain pinned while other parts of the system flow, leading to a plastic 
flow with many defects. One thus need to check again for the stability of the 
moving structure to defects. Fortunately the very existence of channels provide 
a very natural framework to study the effect of such defects: they will lead quite 
naturally to a coupling and a decoupling of the channels ]8~l| , [34 82,8(J. 
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Fig. 34. In a moving glass the motion occurs through elastic channels. The channels 
themselves are rough and meander arbitrarily far from a straight line, but all particles 
follows on these channels like cars on a highway. These channels have been observed in 
decoration experiments, as shown on the left [p3| 



The various phases that naturally emerges in d — 3 are thus the ones shown 
in Fig. ^ (a similar study can be done for d = 2). At large velocity the channels 
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Fig. 35. The various dynamical phases of a moving periodic system. 
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are coupled and the system possesses a perfect topological order (no defects). 
The moving glass system is thus a moving Bragg glass. The structure factor has 
six Bragg peaks (with algebraic powerlaw divergence) showing that the system 
has quasi-long range positional order. If the velocity is lowered a Lindemann 
analysis shows that defects that appear first tend to decouple the channels. This 
means that positional order along the direction of motion is lost, but since the 
channel structure still exist a smectic order is preserved (channels become then 
the elementary objects). As a result the structure factor now has only two peaks. 
This phase is thus a moving smectic (or a moving transverse glass, as first found 
in [jg5| ). It is important to note that in these two phases the channel structure 
is preserved and described by the moving glass equation. Both these phases 
are thus a moving glass. A quite different situation can occur if the velocity is 
lowered further. In that case the channel structure can be destroyed altogether, 
leading to a plastic phase. Note that depending on the amount of disorder in 
the system this may or may not occur, and in d = 3 a purely elastic depinning 
could be possible in principle (in d = 2 the depinning is always plastic). These 
various phases lead to the dynamics phase diagram shown in Fig. pq. These 




Fig. 36. The dynamical phase diagram j82| as a function of the temperature T, the 
magnetic field (far from H c \) and the applied force F. Both the Bragg glass and the 
Moving Bragg glass have perfect topological order. 



various phases as well as the dynamical phase diagram have been confirmed in 
numerous simulations (see e.g. j85| , |8^ , ^ , ^8| ) . 

Let us now turn to the experimental analysis of the dynamics. Early experi- 
ments on the dynamics of the vortex phases near the peak effect (5?|] showed that 
not only does the critical current rise sharply at the transition, but the I — V 
curves also change in character. This is shown schematically in fig. 37. The top 
panel shows the peak effect at a fixed T and varying H . Three distinct types of 
I-V curves are observed in the regions marked I, II and III, shown in the mid- 
dle panel. In the peak region the curves show opposite curvature to that in the 
other regions, the voltage grows convex upwards. The lower panel illustrates this 
behavior through the differential resistance Rd (= dV/df)ior each region. In the 
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Fig. 37. Summary of the variation in the I-V characteristic across the peak effect, 
i..e., ordered to disordered phase transition p^j ]. The upper panel shows the peak effect 
transition and marks the three regimes of flow. The middle and lower panels show the 
I-V curves and the dV/dIT curves for the three regimes. 
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peak region it shows a pronounced maximum above which it rapidly decreases to 
a terminal asymptotic value of the Bardeen-Stephcn flux flow resistance above 
a crossover current. A comparison with simulations shown in Fig. |38| sug- 
gests that the peak signifies a plastic flow region where the vortex matter moves 
incoherently and a coherent flow is recovered at high drives. This behavior con- 




Fig. 38. Simulation data |g9| of the force- velocity curves together with the variation of 
defect density with driving force. Compare the simulation data with the experimental 
data for plastic flow. 
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trasts with that in region I below the peak effect where the pinned vortex matter 
is in the Bragg glass phase. In this case the flux flow resistance is approached 
monotonically without a measurable signature of plastic flow. The plastic flow 
regime is also accompanied with pronounced fingerprint effect: a repeatable set 
of features in Rd as the current is ramped up and down, signifying a repeatable 
sequence of chunks of pinned vortex assembly depinning and joining the main 
flow. This clearly establishes the defective nature of the moving vortex matter in 
this regime. The contrast with this behavior in region I is thought to represent 
a depinning of the Bragg glass phase directly into a moving Bragg glass phase. 

A striking result is seen in the behavior of the flux flow noise j30| , summarized 
in Fig. [39]. The noise is larger in the peak regime by orders of magnitude (seen in 
the lower panel) and is restricted to current values between the depinning current 
and the crossover current (shown in the upper panel). The noise is, therefore, 
associated with an incoherent flow of a defective moving phase and a coexistence 
of moving and pinned vortex phases as was seen in detailed studies |9l|j9^| . 
The qualitative behavior observed in the experiment is consistent with many 
simulations of the flux flow noise characteristics. However, recent work suggests 
that an edge contamination mechanism in the peak regime is responsible for 
triggering much of the defective plastic flow in the bulk and is a subject under 
active investigation ]9^ , p4| . There have also been reports of narrow band noise 
[ f)5| , i.e., noise peaks at the so-called washboard frequency. The observation of 
such noise would be compatible with the moving Bragg glass phase. 

More experimental evidence in favor of the moving glass is also found in 
recent decoration and STM studies of the moving vortex assembly Q,^,^7) . 



5.3 Metastability and history dependence 

A current focus on experimental vortex phase studies is a renewed interest in 
history effects that have long been known to occur in vortex matter. In order 
to understand the equilibrium phase behavior of the system, we need to ascer- 
tain that we have indeed reached equilibrium in an experiment. However, most 
experiments in low Tc systems and at low temperatures in high Tc systems as 
well, show a pronounced dependence of the vortex correlations on the thermo- 
magnetic history of th e sy stem. A striking example is shown in Fig. ^ where 
the magnetic response and transport critical current |l00|] are measured 

for field cooled (FC) and zero field cooled(ZFC) cases. In the latter case one sees 
a pronounced peak effect but not in the former. In other words, the FC state 
yields a highly disordered vortex glass phase and the latter yields an ordered 
Bragg glass phase. The question then is : which is the stable and equilibrium 
state of the system and how does one find it ? A similar question arose in the 



interpretation of neutron experiments [ 101 1 where it was suggested [102] that 



the observed broadening of the lines that disappeared after a cycling above the 
critical current, was due to the presence of out of equilibrium dislocations (on 
the top of the equilibrium Bragg glass phase). 

One possible resolution of the problem comes from a "shaking experiment" 
where the FC state is subjected to a large oscillatory magnetic field and the 
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Fig. 39. Noise characteristics in the plastic flow regime ^(J. The upper panel shows 
the current dependence of the noise which turns on at the critical current and depletes 
at large current when a moving ordered phase is recovered. The middle panels show 
the field dependence of the noise; it is restricted only to the peak regime. The lowest 
panel show a measure of the non-gaussianity of the noise which is maximum (i.e., few 
noise sources) at the onset of motion supporting the moving "chunk" scenario in the 
plastic flow regime. 
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Fig. 40. History dependence of t he c ritical current seen with magnetic susceptibil- 
ity(left) [Q, and transport (right) [[lOC|] , For a very clean sample on the top panel on 
the left, no history effect is seen unlike the bottom panel for a dirty sample. The peak 
of the peak effect marks the onset of an equilibrium disordered phase. 



system evolves to the ZFC state [ 103 1 . Once in the ZFC ordered state, below 
the peak effect, the system cannot be brought to the FC disordered state re- 
gardless of any external perturbation. It is thus reasonable to conclude that the 
ordered (Bragg glass) state is the equilibrium state below the peak effect. The 
FC disordered (vortex glass) phase is simply supercooled from the liquid phase 
above. When pinning sets in, the system fails to explore the phase space due 
to the pinning barrier and stays frozen into glassy, disordered phase. Shaking 
with an ac field then provides an annealing mechanism to bring the system in 
the true ground state which is the Bragg glass. On the other hand, shaking fails 
to produce an ordered state above the peak and thus the disordered phase is 
indeed the ground state there. Curiously then the ZFC state is formed by vor- 
tices entering the system at high velocity, thereby ignoring pinning and forming 
a moving Bragg glass phase from which the pinned Bragg glass phase evolves 
easily. Recent aging experiments 104 have provided compelling evidence in sup- 
port of these conclusions. Further evidence of a thermodynamic nature of the 
transition is obtained also from magnetization anomalies from annealed vortex 
states [ 105 1 . 



5.4 Edge effects 

Yet another phenomenon has begun to be explored in experimental studies 
[ p3|j9^| . This relates to the observation that edges of a sample provide nucle- 
ation centers for the disordered phase. The net results are (1) the order-disorder 
phase transition becomes spatially non-uniform and leads to phase coexistence 
that marks the width of the peak effect region and (2) an external driving cur- 
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rent flows non-uniformly in the system. The contamination of the ordered phase 
by the disordered phase from the edge and the subsequent annealing back to 
the ordered phase at larger drives occur in a non-uniform manner leading to a 
variety of unusual time and frequency dependent phenomena observed earlier. 
Differentiating the effects of these processes from the bulk response of the sys- 
tem, usually assumed in interpreting data and in simulations, remains a subject 
of current study. 



5.5 Transverse critical force 

One of the unexpected consequences of the presence of channels in the moving 
glass phase, is the existence of transverse pinning |j3l| . Let us examine what 
would happen if one tried to push the moving system sideways by applying a 
force F y in a direction perpendicular to motion. This is depicted in Fig. fO] The 




naive answer would be that the system is submitted to a total force F = F x + F y 
and because the modulus of this force is larger than the threshold (F x > F c 
since the system is already sliding), the system will slide along the total force. 
This means that there is a linear response v y to the applied transverse force 
F y . This is what would occur if the moving system was a crystal. However in 
the moving glass although the particles themselves move along the channels, 
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the channels are submitted to a static disorder and thus pinned. It means that 
if one applies a transverse force, the channels will have to pass barriers before 
they can move. As a results the system is transversally pinned even though it 
is moving along the x direction. The existence of this transverse critical force 
is a hallmark of the moving glass. It is a fraction of the longitudinal transverse 
force and decreases as the longitudinal velocity increases as shown on Fig. |d|. 
The existen ce of suc h transverse force has been confirmed in many simulations 
(see e.g. (85,106,107 ). Its experimental observation for vortex systems is still an 



experimental challenge. 



6 Conclusions and perspectives 

It is thus clear from the body of work presented in these notes that the field 
of vortex matter has considerably matured in the last ten years or so. Many 
important physical phenomena have been unravelled, and a coherent picture 
starts to emerge. As far as the statics of the problem is concerned we start 
now to have a good handle of the problem. The Bragg glass description has 
allowed to build a coherent picture both of the phase diagram and of most of 
the previously poorly understood striking phenomena such as imaging, neutrons 
and peak effect or more generally transport phenomena. Important issues such 
as the contamination by the edge and the necessity to untangle these effects 
to get the true thermodynamics properties of the system are now understood, 
allowing to get reliable data. Although, understanding the dynamics is clearly 
much more complicated both theoretically and experimentally, here also many 
progress have been made. The vortex matter has allowed to introduce and check 
many crucial concepts such as the one of creep motion. For the dynamics also it 
is now understood that periodicity is crucial and that strong effects of disorder 
can persist even when a lattice is fast moving. 

Of course the progress realized make only more apparent the many exciting 
issues not yet understood and that are waiting to be solved. For the statics the 
nature of the high field phase above the melting field of the Bragg glass is still to 
be understood. Is this phase a distinct thermodynamic phase from the liquid ? 
is it a true glass in the dynamical sense ? Many difficult questions that will need 
the understanding of a system in which disorder and defects (dislocations, etc.) 
play a crucial role. Similar question occur for the dynamics whenever a plastic 
phase occurs. There is thus no doubt that understanding the role of defects and 
disorder is now one of the major challenge of the field. The glassy nature of the 
various phases also certainly needs more detailed investigations, in particular the 
"goodies" such as aging that have been explored in details for systems such as 
spin glasses would certainly prove useful to investigate. Finally, both theorists 
and experimentalists have mostly focussed for the moment on the steady state 
dynamics. Clearly the issue of noise, out of equilibrium dynamics, and history 
dependence are challenging and crucial problems. 

In addition, the material covered in these notes is only a part of the many in- 
teresting phenomena related to the physics of vortex lattices, and more generally 
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disordered elastic systems. Because of space limitation, it would be impossible 
to cover here all these interesting aspects but we would like to at least mention 
few of them. 

In addition to point like impurities there is much interest, both theoretical 
and practical, to introduce artificial disorder. The most popular is the one pro- 



duced by heavy ion irradiation [108 10E 



defects and to the Bose Glass phase [111]. Other types of disorder such as splay 



11C |, leading to the so called columnar 



[ p_12| or regular pinning arrays have also been explored. We refer the reader to 
the above mentioned reviews for more details on these issues. 

Quite interestingly, the description of the vortices in term of elastic objects 
allows to make contact with a large body of other physical systems who share 
in fact the same effective physics. This ranges from classical systems such as 
magnetic domain walls, wetting contact lines, colloids, magnetic bubbles, liquid 
crystals or quantum systems: charge density waves, Wigner crystals of electrons, 
Luttinger liquids. These systems share the same basic physics of the competition 
between elastic forces that would like to form a nice lattice or a flat interface 
and disorder. This makes the question of exploring the connections with the 
concepts useful for the vortex lattices (similarities/differences) particularly in- 
teresting and fruitful. For example, m agne tic domain walls a re id eal systems to 
quantitativ ely c heck the creep formula |l!3| , Wigner crystals [ 114 1 or charge den- 
sity waves [115 1 have been fields to test for the existence of a transverse critical 
force. The question of the observation of the noise in these various systems is 
also a very puzzling question. Of course, these systems have also their own par- 
ticularities and present their own challenging problems. The interested reader is 
again directed to the review papers for the classical prob lems and a short review 
along those lines for quantum systems can be found in [116|. 

There is thus no doubt that the field is now growing and permeating many 
branches of condensed matter physics, providing a unique laboratory to con- 
tinue developing unifying concepts such as the ones born and grown for vortices 
since more than half of a century. To what new developments this extraordinary 
richness of physical situations will lead, only time will tell. 
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